Reversible limit of processes of heat transfer 
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We study a process of heat transfer between a body of heat capacity C{T) and a sequence of 
A*' heat reservoirs, with temperatures equally spaced between an initial temperature To and a final 
temperature Tjv. The body and the heat reservoirs are isolated from the rest of the universe, and 
the body is brought in thermal contact successively with reservoirs of increasing temperature. We 
determine the change of entropy of the composite thermodynamic system in the total process in 
which the temperature of the body changes from To to Tjv. We find that for large values of A'^ the 
total change of entropy of the composite process is proportional to (Tjv — To)/N, but eventually a 
non-monotonic behavior is found at small values of A'^. 

PACS numbers: 



I. INTRODUCTION 

One of the original formulations of the second law of 
thermodynamics, due to Clausius, states that "no pro- 
cess is possible whose only result is the transfer of heat 
from a body of lower temperature to a body of higher 
temperature". As discussed by Fermi in his book [l|, 
since as this statement was made the concept of temper- 
ature was empirical, to find out which of two bodies has 
the higher temperature one had to put them in thermal 
contact and find out in which sense the heat flows, it will 
flow spontaneously from the body of higher temperature 
to the other. So, Clausius's statement could be rephrased 
into "If heat flows spontaneously from body A to body 
B, no process is possible whose only result is the transfer 
of heat from body B to body A" . To be more precise, 
let us isolate the two bodies from the rest of the uni- 
verse, to assure that they will not exchange heat or other 
forms of energy with the environment, this is implicit in 
the term 'only result'. We may notice in this statement 
that already the idea of irreversibility is present: if heat 
flows spontaneously from A to B, it will not flow spon- 
taneously from B to A. In other words, the process of 
transfer of heat between two bodies which are isolated 
from the environment is irreversible. 

The second law off thermodynamics says that not all 
processes which conserve energy, as required by the first 
law, occur in nature. In many cases (if not all), if a given 
process in an isolated system occurs, the inverse process 
does not. This idea was further formalized by the intro- 
duction of the concept of entropy, so that processes in 
isolated systems which lead to a decrease of the entropy 
are not allowed. In particular, in the process of heat 
transfer mentioned above, the entropy of the composite 
isolated system increases if the heat flows from the body 
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of higher temperature to the one of lower temperature, 
but would decrease in the inverse process, so that only the 
further is allowed. It is interesting to recall that in his 
paper of 1885, Clausius summarized the two first laws 
of thermodynamics in just two sentences: "The energy 
of the universe is constant. The entropy of the universe 
tends to a maximum" 0. We may notice that the uni- 
verse is the only composite thermodynamic system which 
is isolated by definition. 

The entropy in equilibrium thermodynamics is a state 
function, so that it is defined for equilibrium states of 
the system. When two bodies of different temperatures 
are brought into thermal contact, the composite system 
in general will go through a sequence of non-equilibrium 
states, but if we wait long enough, it will finally reach 
a new equilibrium state, where both bodies of the com- 
posite have the same temperature, and in this state the 
entropy is larger than it was in the initial state. If we now 
imagine a process which happens at a very slow rate in 
time when compared to the relaxation time of the system 
Q , we may reach a situation close to the one in which all 
intermediate states of the system are equilibrium states. 
In this limit the process is called quasi-static, and there- 
fore the whole process may be represented as a path in 
the space of variables which define the equilibrium states 
of the system. In the particular case in which the entropy 
of the final state of the quasi-static process is equal to the 
entropy of the initial state, and therefore to the entropies 
of all intermediate states, this process will be reversible. 

The process we will study in some detail here is the 
transfer of heat in isolated systems composed by a heat 
reservoir (constant temperature) and a body whose heat 
capacity is described by a function C(T). The initial 
temperature of the body is Tq and its final temperature 
is Tjv. The change of the temperature of the body is 
accomplished by putting it in thermal contact with a 
sequence of N heat reservoirs whose temperatures are 
equally spaced between Tq and Tpf. For a finite num- 
ber of reservoirs, each of the sub-processes consists of 
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heat transfer between two bodies at different tempera- 
tures and is therefore irreversible, leading to an increase 
of the entropy of the composite isolated system. In the 
limit N oo, however, each sub-process will contribute 
with a negligible increase of entropy and, what may not 
be obvious initially, the total variation of entropy van- 
ishes as well, so that this is a concrete example of a limit 
which leads to a reversible, quasi-static process of heat 
transfer. 

This problem, of considerable pedagogical interest, 
has been studied before in the literature. For a con- 
stant heat capacity, such as is found for ideal gases, 
Calkin and Kiang [4] have shown for a cyclic process 
(To —7' T/v — !> To), that the total change of entropy van- 
ishes as 1/N in the limit N ^ oo because, although 
the number of sub-processes increases, the increase of 
entropy in each of them is proportional to in this 

limit. A similar situation is proposed in the exercise 4.4-6 
in Callen's book Q. The same process was studied also 
by Thomsen and Bers fs'], again for the case of constant 
heat capacity and discussing in detail that it is always 
possible, given any positive value e, to choose a value of 
N for which the increase of entropy is smaller than e. 
This argument is then repeated for the inverse process, 
thus characterizing the process in the limit iV — !■ cx) as 
reversible. 

Here we generalize the problem in two ways: the heat 
capacity of the system is no longer constant and we de- 
velop the increase of entropy as a series in 1/iV, so that we 
find the corrections to the leading term. We also show 
that some symmetries in the result which are valid in 
the asymptotic limit close to reversibility are not longer 
present for finite values of A''. In particular, we show that, 
although in the large N limit the increase of entropy de- 
creases monotonically with N, this may not happen at 
lower numbers of heat reservoirs. Also, the increase of en- 
tropy for the direct process (Tq Tn) may be different 
of the same quantity for the inverse process (T/v Tq). 
The determination of the change of entropy is given in 
section |lT] and applications to particular cases and final 
comments may be found in section Hill 



II. DETERMINATION OF THE CHANGE OF 
ENTROPY 

Let us define the problem in more detail. Our com- 
posite system consists of the body, whose heat capacity 
is C(T) and which will undergo a temperature change 
AT = T/v — To, and a set of N heat reservoirs, at 
equally spaced temperatures, so that the temperature 
of the reservoir j {j = 1,2,3, N) is equal to Tj = 
Tq + j AT/N. The body is brought into thermal contact 
with the heat reservoirs in order of increasing values of j 
until its temperature equals the one of the reservoir and 
therefore is increased bi AT/N. In the first sub-process, 
for example, the body starts with temperature Tq and 
ends with temperature To + AT/N. After the TV'th pro- 



cess, the body will reach the final temperature T/v. This 
sequence of heat transfer processes is illustrated in fig- 
ure [TJ Actually, one could argue that the time required 
to reach equilibrium for each sub-process of heat transfer 
would be infinite, and for this reason, Thomsen and Bers, 
in their discussion of a similar process Q , have chosen the 
temperature of the j'th reservoir to be To -I- (j -I- 1)AT/N, 
thus assuring that the equilibration time of the system 
with each reservoir is finite. Here, for simplicity, we will 
not use this improved version of the process. 



FIG. 1: Illustration of the j'th process of heat transfer, where 
the system s is in thermal contact with the heat reservoir at 
temperature Tj, so that its temperature is changed from Tj-i 
to Tj. The pair system-reservoir at Tj is isolated, as are all 
the other reservoirs. 

If a the body, at temperature T, receives a quan- 
tity of heat dQh, its temperature will change by dT — 
dQh/C{T). The change of entropy of the body in the 
whole sequence of heat transfer processes is: 
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The change of entropy of all reservoirs in the sequence of 
processes will be: 
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where in the last passage, we recall that in each sub- 
process, the body and the reservoir are isolated, so that 
the heat received by the reservoir is given by dQr = 
—dQb = —C{T) dT. We thus may write the total change 
of entropy as: 



AS{N) = ASb{N) + ASr{N) = 
''^CiT),^_J2 



To 



T 



C{T) dT.{2) 



If it is possible to perform both integrations, this expres- 
sion above becomes: 



AS{N) = 5fc(T/v) - SbiTo) - ^ ^[U{T,) - C/(T,_i)], 

(3) 
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where Sb{T) = JC{T)/TdT and U{T) = JC{T)dT. It 
is convenient to rewrite the expression for the change of 
entropy as: 



N 



n,{T,N)4>{T) dT, 



(4) 



where Kj{T,N) ^ 



1 - T/Tj and (j>{T) = C{T)/T. It 
is now easy to see that, since stabihty imphes that 
4>{T) > 0, we must have AS* > 0, as required by the 
second law of thermodynamics. Let us show this in some 
detaiL Suppose that Tj\j > Tq, so that the temperature 
of the body increases in the process. It is then clear that 
Kj{T,N) > and as a consequence the entropy of the 
composite system increases. The same happens in the 
inverse process, where the temperature of the body de- 
creases from T/v to Tq. In this case the change if entropy 
is: 



AS'iN) 
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«j_i(r,iV)</>(T)dT, (5) 



and since now Kj_i(T, N) < in the range of integration, 
again the entropy increases. It may be worth noticing 
that the difference of the entropy increases between the 
processes with raising and lowering of the temperature of 
the body SS{N) = AS'{N) - AS{N) is: 
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0(r)dT, (6) 



and the sign of this expression is not defined in general. 

Another general property of the entropy increase is its 
dependence of the number of reservoirs N. As we will 
see below, it vanishes in the limit — > oo of a quasi- 
static process. The question we will address is if the 
change of the entropy decreases monotonically with N. 
For convenience, we will define the function: 

N 

K{T,N) = J2e{T~ T,_i)e(r, - t) k,{t, n), (7) 

i=i 

where 8(r) is the step function. The function k{T,N) 
is defined in the whole temperature range [^o^ ^a^] and 
has a sawtooth pattern, with N maxima of decreasing 
values. The expression for the entropy increase Eq. ([4]) 
may then be cast into the following form: 



AS{N)^ / ^iT,N)m 

J To 



dT. 



(8) 



The difference of entropy increases for two values of 
iN2 > Ni) is: 



y 



FIG. 2: The function n{T, 5) — k{T, 6) for initial temperature 
To — 10 and final temperature Tn = 20. 



It is straightforward to see that if N2 is a multiple of A^i 
the difference of k functions is non-negative, and there- 
fore AS{Ni) — AS{nNi) > for any integer n. For exam- 
ple, we have AS'(l) > AS{2). This, however, is no longer 
true in general if the ratio N2/N1 is not an integer. As an 
example, in figure [2] we show the difference of k functions 
for A^i = 5 and A^2 = 6, and we notice that this differ- 
ence assumes negative values in part of the temperature 
domain. We thus reach the conclusion that, depending 
on the function 0(T), the entropy increase for N2 heat 
reservoirs may be larger than the one for A^i reservoirs, 
if 1 < N2/N1 < 2. 

We now proceed obtaining the general asymptotic be- 
havior of AS{N) ioT N ^ 1. Expanding the last integral 
in the expression for the increase of entropy, Eq. ([2|), 
which corresponds to the amount of heat Qj exchanged 
by the body with reservoir j, in a Taylor series around 
the temperature Tj, we get: 
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where a superscript between parenthesis means a deriva- 
tive of the corresponding function of the order of the su- 
perscript. We may now write the total change of entropy 



AS'(A^) 



C{T) 
T 



dT 



T, 



(11) 



AS{N^) - AS{N2) = / [k(T, A^i) - k(T, iV2)]0(T) dT. 
JTq 

(9) 



We now proceed transforming the sum over the reser- 
voirs into an integral, using Euler-MacLaurin's expansion 
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where i?2fc are the BernouUi numbers {B2 ~ 1/6, i?4 = 
-1/30, Be = 1/42, = -1/30, . . .. See [7] for a hst of 
the numbers). The series in the right hand side is diver- 
gent in many cases, since the Bernoulh numbers increase 
very fast for larger values of k. There are procedures to 
evaluate the rest Rkr„,a,^ l6| and we will discuss this issue 
in the appendix. The sum we will convert into an integral 
is: 



^C('-i)(T,-) 
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To 



and applying the expansion (jl2[) we find: 
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Now we may substitute this last expression into Eq. 
(fTTjl and collect the terms in the change of entropy of 
the body in the sequence of processes in powers of 1/A^. 
We notice, as expected, than the term of order cancels, 
so that, in the limit N ^ 00 the whole process becomes 
reversible. In the appendix, we proceed calculating all 
coefficients ai of the expansion 
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but here we we will limit ourselves to the leading term, 
which is: 
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(16) 

Integrating by parts, we finally obtain the leading term 
of the expansion: 



AS'(iV) 



Tn — To 
2N 



Tn 



CjT) 

rp2 



dT, 



(17) 



valid in the limit iV ^ 1. 

Since G(T) > 0, we notice that the leading term van- 
ishes only in the trivial case where G (T) = in the whole 
range of temperatures and there is no heat transfer, so 
that the asymptotic behavior as 1/iV in the increase of 
entropy, as was found in the particular case of an ideal gas 
U, is universal. We also remark that in the asymptotic 
regime the entropy increase is invariant if the tempera- 
tures are interchanged and also monotonically decreasing 
with N , properties which are not true in general for small 
values of the number of reservoirs. 



III. DISCUSSIONS AND CONCLUSION 

We will now apply the results of the preceding section 
to some particular cases of physical interest. We start 
with the ideal gas, which has a constant heat capacity 
G(T) = Co- Applying Eq. ^ to this particular case, the 
increase of the entropy will be: 
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where / = Tn/To- The leading term in this case is: 

In figure [3] we show some curves of the entropy increase 
as a function of 1/A^, for different values of the ratio 
of temperatures /, as well as the dashed straight lines 
which represent the asymptotic behavior. Notice that, 
as expected, the same asymptotic behavior is found for 
/ = 2 and / — 1/2. If we consider two processes, one 
with ratio /i > 1 and the other with /2 = I//1 < 1, that 
is, related by interchange of the temperatures, both will 
have the same asymptotic behavior, as may be seen in 
the particular cases depicted in the figure. Also, the in- 
crease of entropy in the second process, where the body 
is cooled, is always larger than the one in the first pro- 
cess. Another way to state the same property is that, 
as can be seen in the curves, the increase of entropy is 
a concave function of 1/iV for processes where the body 
is heated and it is a convex function when the temper- 
ature of the body decreases. In the asymptotic regime 
these properties may be obtained the next term of the 
Euler-MacLaurin expansion, as will be discussed in the 
appendix. Also, the increase of entropy for this particular 
case is a monotonically decreasing function of N. 

In most cases the variation of entropy decreases mono- 
tonically with N, a property which seems also intuitive, 
since by decreasing the temperature steps in a certain 
sense the composite process gets closer to a reversible 
process. For the ideal gas this is always the case. We have 
seen above, however, that this may not be generally true. 
If we look at the Eq. © for the difference of entropy in- 
creases for different values of A^, a non-monotonic result 
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FIG. 3: Increase of entropy for a system with constant heat 
capacity Co (ideal gas) as a function of the inverse number of 
reservoirs 1/A''. The four curves, in the upward order, corre- 
spond to ratios of temperatures / — Tm /To — 2, 1/2, 5 and 
1/5. The dashed straight lines are the asymptotic behavior 
for small 1/N. 

is possible if the function 0(T) has one or more peaks at 
the temperatures where k{T,Ni) — k{T,N2) is negative. 
One system which could possibly show a non-monotonic 
behavior is a two-state system [3|, since its heat capac- 
ity displays a peak, sometimes called "Schottky hump" . 
However, we found that this is not the case, the reason 
is that the peak is too wide to produce a non-monotonic 
behavior. 

We turn our attention to a situation where the width 
of the peak in the heat capacity may be changed by ad- 
justing a parameter. The simple expression we will use 
for the heat capacity is: 

C{T) = C[t{2 - r)]"\ (20) 

where r = T/Tc, so that Tc is the temperature where 
the maximum of the heat capacity is located, and C is 
a constant with the dimension of entropy. Although this 
heat capacity does not correspond to a particular physi- 
cal system, it is convenient for analytic calculations and 
consistent with the third law of thermodynamics. Of 
course, at least for odd values of the parameter m, it is 
valid only if r is in the range [0,2], so we will restrict 
ourselves to this range of temperatures. As the param- 
eter m grows, the peak in C{T) narrows, and our main 
interest in this system is that, for sufficiently large values 
of m, non-monotonic behavior is found in AS{N). 

In figure 0] we show, in the same graph, the function 
k(t, 2) — k(t, 3) and the function 0(t) = C(t)/(Ct), 
given by Eq. (j20p with m — 16. As is apparent in the 
graph, with the choices of initial (tq = 0.286) and final 
(r^v = 2.0) temperatures we made, the maximum of the 
peak in C (T) is close to the center of the negative peak 
in k(t, 2) - K(r, 3), so that the difference AS (2) - AS{3) 
is minimized, since this difference is given by Eq. ([9]) . In 
tableUwe list the resuhs for [A5(3) - AS{2)]/C for dif- 
ferent values of m, with the choices above for tq and r^r, 
and it is clear that a non-monotonic behavior is found 



FIG. 4: k(t, 2) — k(t, 3) (dashed line) and the function <^(r) = 
C{t) / {Ct) = [r(2 — t)]^^ /t. The initial temperature is ro = 
0.286 and the final temperature is rjv = 2.0 



m 


[AS'(3) - AS(2)]/C 


14 


-3.158309 X 10"^ 


15 


-3.427008 X IQ-* 


16 


2.185200 X IQ-^ 


17 


4.463576 x 10"^ 



TABLE I: [AS'(3) - AS'(2)]/C for difi'erent values of the pa- 
rameter m for a system with heat capacity given by Eq. (I20|l . 
The temperature range of the process is ro = 0.286, tm = 2.0. 

for m > 16. Finally, we found some evidence that a 
local maximum in the entropy could also be found at 
= 4, but if this actually happens for this particular 
form of heat capacity, it will be at rather large values of 
TO, which lead to numerical problems. We remark that, 
in order to reduce numerical errors, in all examples we 
discuss the integrations involved in the calculation of the 
entropy increase Eq. ^ are performed analytically, as is 
done in Eq. 

As discussed before, we expect, qualitatively, the non- 
monotonic behavior to be enhanced if the peak in the 
heat capacity is more pronounced. This leads us, in this 
final example, to a system which undergoes a continuous 
phase transition, since in this case a singularity at the 
critical temperature Tc is found in the heat capacity, of 
the form C{T) « |1 - T/Td"", described by the critical 
exponent a [8|. To the singular behavior of the heat 
capacity, usually regular contributions should be added, 
so that we will use the simple expression: 

C(T) = Cr|l-r|^", (21) 

where we define the reduced dimensionless temperature 
T T/Tc and which is compatible with the third law 
of thermodynamics. It should be mentioned that since 
the exchange of heat when the temperature of the body 
crosses the critical temperature has to be finite, we should 
restrict the critical exponent a to values smaller than 1. 
It is easy to perform the necessary integrations in this 
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FIG. 5: Values of AS{N) as function of 1/7V for a body with 
a continuous phase transition (Eq. (|21|l '). The upper curve is 
for Q = 0.7 and the lower one for a — 0.5. The pairs of data 
at the end of the thick lines are the first ones, in the order of 
increasing values of A'^, where a increase of AS{N) is found. 
These results are for tq = 0.2 and tjv = 2.2. 



case, so that the entropy increase will be given by Eq. 
dS]) with: 



c 



(l-r) 



1-a ' ^< 1' 
T> 1, 



(22) 



and 



(l-a)(2-a) 



(1 - a)r], r < 1 



(23) 



r > 1. 



We would expect that the non-monotonic behavior of 
the entropy should be enhanced as the value of the critical 
exponent a increases. This is actually true, but at vari- 
ance with what is seen in the preceding example of a peak 
in the heat capacity, the local maxima of the entropy in- 
crease, for lower values of a, appear at larger values of 
N . In figure [S] we see results of the increase of entropy 
for two values of the exponent, a = 0.7 and a = 0.5. It 
is apparent that, for a = 0.7 we already have an increase 
of the change of entropy between N\ — 2 and A^2 = 3, 
the segment with a thicker line in the curve. Many other 
pairs of successive results with the same property are 
found for larger values of TV. For a = 0.5, the first pair 
where the change of entropy increases is seen at A^i — 7, 
N2 = 8. 

For a given a value of the critical exponent a, we may 
find A^i, the lowest value of N such that AS{N2) > 
AS{Ni), where A^2 = ^1 + 1- In figure [BJ the exponent 
a is plotted as a function of 2/(A^i -I- A^2), and we notice 
that, at least for large values of a, that the function shows 
steps. The results shown are calculated for tq = 0.2 and 
tn ~ 2.2. Thus, for example, if a > 0.555, we notice 
that the first increase of AS{N) is seen between Ni = 2 
and = 3, while if a is in the range [0.555,0.454] this 



FIG. 6: Values of the critical exponent a for which the first 
increase of AS{N) if found between Ni and = A^i + 1, 
as function of 2/{Ni + N2). The data are for ro — 0.2 and 
rjv = 2.2. 



increase is seen at iVi = 7, N2 — 8, if a is further de- 
creased the increase shifts to iVi — 12, N2 — 13, and so 
on. This pattern, although qualitatively remaining the 
same, changes quantitatively for a different temperature 
range. We notice an interesting behavior for small val- 
ues of a, which leads to the question if, for any positive 
a an increase if AS{N) is found for a finite value of N. 
Due to round-off errors, we could not address this ques- 
tion numerically. However, since the coefficients ai in the 
Euler-MacLaurin expansion Eq. (I15p diverge in this case 
for i > 1 and any positive a, the answer to this question 
should be positive. We discuss this point somewhat more 
in the appendix. 

Finally, we notice that, in opposition to what was 
found in the first two examples we studied, when the 
body undergoes a continuous phase transition it is pos- 
sible that AS{N) > AS'{N), that is, the increase of en- 
tropy in the heating process is larger than the one in the 
cooling process, a possibility which was anticipated in the 
discussion after Eq. ^ . In figure [7] we present results 
for the difference in entropy increases AS'{N) — AS{N), 
for a — 0.7 and a = 0.5, as functions of The tem- 

perature interval is tq — 0.2, tn = 2.2 Although in most 
cases this difference is positive, for a = 0.7 already for 
n = 13 we get a negative result, and a set of values of N 
above this one have the same property. For a = 0.5 the 
effect is smaller, and the first occurrence of a negative 
result is for = 88. Again it would be interesting to 
study more details of these results at large values of TV, 
but round-off errors prevent this to be done numerically. 

A final remark on this example is that, if r <i; 1, 
C(r) w Cr, so that it will show a behavior similar to 
a gas of non-interacting fermions at low temperature. In 
this case, if the initial temperature tq approaches 0, the 
coefficient cri diverges as a logarithm. In other words, in 
a heating process starting at tq = 0, the tangent to the 
curve AS* x 1/iV at the origin is vertical. This curve is 
concave, since (T2 is negative. The even more unphysical 
cooling process with vanishing final temperature shows a 
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FIG. 7: Difference in the increases of entropy in the coohng 
and the heating processes AS'{N) — AS{N), for to = 0.2 and 
Tjv = 2.2. The curve with larger oscillations is for a = 0.7 
and the other one for a = 0.5. 



very uncommon behavior: AS{N) diverges for any finite 
N and vanishes as — oo. 

In conclusion, we notice that in general the process of 
heat transfer between a body and a sequence of heat 
reservoirs, is a very rich physical situation for under- 
standing of the second law of thermodynamics and the 
concept of a reversible process. In particular, if the body 
undergoes a continuous phase transition in the composite 
process, some curious phenomena appear , which at first 
seem to defy common sense, but of course none of them 
violates thermodynamics. 



Appendix: Euler-MacLaurin expansion 

Here we will develop the terms of the Euler-MacLaurin 
expansion in some detail. If we substitute Eq. (jl4[) into 
Eq. ((TT|) and collect the terms in the change of entropy of 
the body in the sequence of processes in powers of AT/N. 
We notice, as expected, than the term of order cancels, 
so that, in the limit N oo the whole process becomes 
reversible. The coefficients we obtain for the change of 
entropy Eq. ([15]) are: 



-1 /•'^" CW(T) 



To 



[ (* + 1)JTo 

1 ( c^'-^Htn) _ c('-i)(ro) 

^ (i + 1 - 2k)\(2k)l^ 

k—1 

Tn ) 



T 



dT+ 



(A.l) 



The two first coefficients of this series, after integrating 
by parts, are:, are: 



TWO r-"£(T) 
2 L r2 



(A.2) 
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{Tn - To) 
12 



C{Tn) CiTo) 



To 



cm 

T3 



dT~ 



^N 



rp2 

^0 



(A.3) 



We notice that ai > 0, and vanishes only if the tempera- 
ture interval is zero or if C(T) = in the whole interval. 
In other words, if there is a non-zero heat exchange be- 
tween the body and the reservoirs, ai is positive. If this 
coefficient could be negative, the second law would be vi- 
olated. It is also invariant with respect to a permutation 
of the temperatures, as was already noted in the exam- 
ples discussed above. The sign of the second coefficient 
tT2 is not fixed, and, as happens with all coefficients of 
even order, it switches under permutation of the temper- 
atures. The odd order coefficients do not change their 
sign under this permutation. 

In the particular case of an ideal gas, where C(T) = 
Co, the coefficients are: 



<7l 



0'2fe 



-Co 



T, 



N 



To 



Co 



2 

B2k {Tn 



1 

TV 
To^'^ 



(A.4) 



2k 



1 



^N -'O 



where the last expression is valid for k=l,2,3,. . . . We no- 
tice that, with the exception of the dominant term, only 
terms with even powers of 1/A^ appear in the expansion 
for this particular case. For the heating process Tn > To, 
(72 < 0, so that the AS" is a concave function of 1/iV in 
the limit N ^ oo. For the cooling process Tn < To, the 
function is convex. These properties can be verified in 
the numerical results presented in figure [3] 

Very often the expansion in the Euler-MacLaurin sum- 
mation (IT^ does not converge. We find that even for 

= 1 the truncated asymptotic values are close to the 
exact ones at rather low orders. However, truncating at 
higher orders does not necessarily lead to results closer 
to the exact ones. For example, for = 1 the result 
closest to the exact one is found truncating the series at 
kmax — 6. We notice that higher orders of truncation 
lead to wrong results, as may be seen in figure [51 where 
the increase of entropy calculated using the expansion 
Eq. ([m, for TV = 1 and Tn/To = 2, is shown for sev- 
eral orders of truncation. The results oscillate around 
the exact value at successive orders. The best result is 
obtained for imax = 6 and if we truncate the expansion 
at larger orders we observe that they start deviating from 
the exact value. 

In the table [TTl we present more results of such calcu- 
lations. We notice that quite accurate estimates of the 
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FIG. 8: Increase of entropy for a system with constant heat 
capacity Co as a function the order of truncation of the asymp- 
totic expansion imax- The ratio of temperatures Tn/Tq is 
equal to 2 and the number of reservoirs is TV = 1. The exact 
value corresponds to the dashed line. 



entropy change may be obtained if the series is truncated 
at the most favorable order. The question of determin- 
ing the order which minimizes the rest R in the Euler- 
MacLaurin expansion Eq. (fT2|) is rather technical and we 
refer to the specialized literature for details, but in gen- 
eral it may be assured that, if f{x) in Eq. ([12]) has always 
the same sign and if the function and all its derivatives 
tend monotonically to as a; — oo, which are valid for 
f{N) = AS{N), then the rest i?fc„„^ is of the same order 
and has the same sign of the first neglected term [Q] . 

Finally, it is clear from Eq. (jA.ip for the expansion 
coefficients, that for the case of a singular heat capacity 
at a critical temperature Tc {C{T) = Ct(1 — r)", with 
T = T/Tc), all but the first coefficients diverge if the 
temperature interval includes the critical temperature. 
This fact may account for the rather peculiar behavior 
found numerically in the large N limit. 
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TABLE II; Body with constant heat capacity {C{T) = Co). 
Values of the relative difference between the increase of en- 
tropy calculated directly (AS) and the one which follows trun- 
cating the series Eq. (|15|l (AS^^^^ ). Several values for Tjv /To 
and A'^ are considered, and the optimal order of truncation 
imax is given. 



